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ABSTRACT 

Renormalization group flow equations for scalar A$ 4 are generated using three classes 
of smooth smearing functions. Numerical results for the critical exponent v in three di- 
mensions are calculated by means of a truncated series expansion of the blocked potential. 
We demonstrate how the convergence of v as a function of the order of truncation can be 
improved through a fine tuning of the smoothness of the smearing functions. 
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I. INTRODUCTION 



Renormalization group (RG) methods provide a powerful tool for investigating non- 
perturbative physical phenomena [1]. Issues such as QCD under extreme conditions, 
formation of a quark-gluon plasma in relativistic heavy-ion collisions [2], or critical phe- 
nomena in condensed matter systems all cannot be treated using standard perturbation 
theory due to the presence of infrared (IR) singularities. Through the continuous elimina- 
tion of degrees of freedom, RG techniques systematically resum the perturbative series and 
therefore can provide information about nonperturbative effects in these systems. However, 
the power of RG relies on the existence of efficient analytic and computational methods 
since the full RG flow equations cannot be solved exactly. The goal of this paper is to 
show how to optimize the way the degrees of freedom are eliminated in order to improve 
results obtained using approximate RG flow equations. 

For a field theoretical system a continuous RG transformation can be realized by 
introducing a smearing function pk(x) which governs the coarse-graining procedure [3] 
[4] . The scale k acts as an effective IR cutoff that separates the low- and high-momentum 
modes. Using this smearing function an averaged blocked field can be defined as 



<Pk(x) = / p k {x - y)<f>(y), 
from which one obtains the effective Legendre blocked action: 



(1.1) 



-s k [*(*)] 



(1.2) 



In this manner we achieve a smooth interpolation between the bare action Sa[3>] defined 
at the ultraviolet (UV) cutoff, A, and the quantum effective action iSfc[<&], which generates 
those one-particle-irreducible graphs whose internal momenta extend between k and A [5] . 
The introduction of a smearing function leads to a modification of the bare propagator: 



A(p) 



A fc (p) 



i-Pfc(p) hip) 



p 



(1.3) 



with pkip) + Pk(p) = 1- One can then derive a RG equation for the blocked action which 
remains valid in all orders of the loop expansion by varying Sk infinitesimally with k: 
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Since the pioneering work of Wilson [1], similar RG equations have been derived and 
analyzed. A sharp momentum cutoff was first used by Wegner and Houghton [6], and 
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further discussed in Refs. [7], [8] and [9]. Polchinski [10], on the other hand, employed a 
smooth cutoff. There exists a vast amount of literature devoted to this subject [11]. Such 
functional RG equations, as stated before, have proven too difficult to be solved exactly, 
and further approximations are needed. Any viable scheme must not only retain the 
nonperturbative characteristics of the RG, but also converge sufficiently rapidly without 
inducing further spurious effects. As demonstrated by Morris in his seminal papers [5] [12], 
the most reliable method so far for probing the low-energy effective theory is the derivative 
expansion: 

= /J^(^) 2 + ^($) + o(a 4 )}, (i.5) 

where and £/&($) are, respectively, the wavefunction renormalization and the blocked 

potential. At leading order in the derivative expansion Zfc(<&) is taken to be unity and the 
low-energy effective action is described solely in terms of £/&(<&) [13]. This local potential 
approximation (LPA) results in the following flow equation: 

k dum i f( k 9pk(p)\ Km 



) (16) 



dk 2j p \ dk J p 2 + pk(p)U'^) 

Although the above expression can be derived exactly [14], it may also be obtained by a 
differentiation of the perturbative one-loop result 

Ui 1 \$) = lj^{l + p k (p)^^} (1.7) 

with respect to ln(k) followed by a substitution V"(&) — > U^{<&) on the right-hand-side. 

Clearly, the functional form of the RG equation for Uk(&) depends on the choice of 
the smearing function. A sharp cutoff, Pk(p) = 0(& —p), provides a well-defined boundary 
between the high- and low-momentum modes and yields a nonlinear partial differential 
equation. On the other hand, for a general smooth cutoff, no clear separation exists 
and the RG flow remains an integro-differential equation. According to the universality 
principle the shape of the smearing function, pk(p), does not influence the physics of finite 
length scales. So long as the effective action contains all marginal or relevant operators 
the resulting RG flow must be scheme independent. At leading order in the derivative 
expansion, with no further approximation involved, physical results are indeed independent 
of the shape of Pk{p), as demonstrated in [5] [15]. However, at next-to-leading order in the 
derivative expansion where the effect of Z^($) is included, the critical properties of the 
system will depend on whether pk(p) is sharp or smooth. In fact, ambiguities arise in the 
former, although self-consistency arguments can be used to circumvent this difficulty [16]. 
The dependence of rj, the anomalous dimension, on the shape of the smearing function 
beyond leading order, has been calculated in [15] [17]. 

The critical exponents and other universal properties of the system may vary when 
they are determined in the vicinity of the fixed point of the truncated equation given in 
Eq. (1.6), or when the theory moves far away from the critical manifold. In fact, even 
if the solution includes all the relevant operators, the neglected irrelevant parameters are 
still evolving at the fixed point of the truncated equation. Since the modification of the 
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irrelevant operators amounts to a change of an overall scale factor of the theory this scale 
factor remains cutoff dependent at the approximate fixed point. In order to minimize its 
variation one has to retain in the approximation the irrelevant operators with the critical 
exponent. 

A frequently used technique is to expand in power series of $ followed by a trun- 

cation at some order, thereby turning the problem into solving a set of coupled nonlinear 
ordinary differential equations. Although the numerical algorithms are simplified by this 
expansion, the convergence of critical exponents calculated in this way is not guaranteed. 
In fact, when employing a sharp smearing function while is expanded in $ about the 

origin, the critical exponent v has been shown to oscillate about its expected value as M, 
the number of terms in the series, is increased [12]. The convergence improves, however, 
if the expansion is made around the /c-dependent minimum of Uk(&) [9]. Unfortunately, 
this approach fails when applied to an O(N) field theory in the symmetry- broken phase 
due to the presence of the massless Goldstone modes which cause IR divergences to persist 
[14] [18]. It is therefore desirable to identify a computational method which does not suf- 
fer from these drawbacks. In addition, for gauge theories, the numerical complication of 
solving the non-truncated flow equations is even more overwhelming because of the prolif- 
eration of degrees of freedom and the complexity of their interactions. For these theories 
a polynomial truncation of the blocked potential seems inevitable. Thus, it is crucial to 
understand how to improve results obtained from truncation schemes. 

In the present work, we examine what happens, in the truncated polynomial expan- 
sion of when a smooth smearing function Pk,a(p) with a being the smoothness 
parameter, is utilized instead of Q(k —p). In particular, we explore how v changes with 
M as the smooth smearing functions approach the sharp limit, and how its convergence is 
influenced by a. Our goal is to obtain a prescription which eliminates or diminishes the 
oscillations mentioned above and ensures a rapid convergence as M is increased. 

The organization of the paper is as follows. In Sec. II we give three examples of 
smooth smearing functions and derive the corresponding RG equations. We show how 
the expected sharp cutoff limit can be recovered by recognizing that Pk(p) is a continuous 
function in the vicinity p ~ k. Finite-temperature RG equations are also derived. In 
Sec. Ill we present the numerical solutions for v associated with the three-dimensional 
Wilson-Fisher fixed point at various levels of polynomial truncation of Uk(&). In Sec. 
IV the source of scheme dependence on M and a is discussed. We propose an optimized 
smooth smearing function which leads to a maximal cancellation of the effects of irrelevant 
operators and provides the fastest convergence of v to a value which is in good agreement 
with the world's best estimate. Sec. V is reserved for summary and discussions. 

II. SMOOTH CUTOFF FUNCTIONS 

We begin by considering a general smooth representation of the smearing function 
Pk,e(p) = ®e(k,p) which approaches Q(k — p) as e — > 0. Due to the singular nature of 
the step function ambiguities can arise when taking the limit e — > 0. More precisely, one 
will encounter terms involving O(0) = 9q. Since this value is not uniquely specified for 
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the smooth parameterizations of the step function, it would seem that the approach to the 
sharp limit is not unique. The naive, or mean, approach in the limit e — > is: 

POO /T, \ POO 

J dp dk G ( Q «(M>p) — J o d P 5(k-p)G(eo,p) = G(e ,k), (2.1) 

where <9fcO e (£;,p) — > — p). This prescription, however, is oversimplified because the 
integrand of 

/ dpG(0 e (fc,p),p)), (2.2) 
Jo 

is not uniformly convergent as e — > 0. The proper treatment requires that e (A;,p) be 
used as an integration variable. Thus by setting 6 e (/c,p) = 6 e (l — p/fc) = t and making a 
change of variables from p to t, we have 



fdp ^^G(e £ (k,p),p) — dp(-l^)G(t,p) = -j^^dt ^ G(f,p(f)). 

t( P -0) ^ 

The sharp cutoff limit is obtained when t(p = 0) = 1, t(p = oo) = 0, p(t) = k and Eq. 
(2.3) simplifies to 

dp d °f k P) G(e £ (k,p),p) = j\t G(t,k). (2.4) 
In terms of t, the RG flow equation in Eq. (1.6) reduces to 



h dU k m = Sd f t(p=oo) df P d (t) U'£{*) 



Ok 2 J t{p=0) p->{t) + {l-t)U>>{<S>) 
S d k d f t(z=oo) z d {t) 



(2.5) 



t( ,=o) z*{t)+{i-t)uzm 



where S d = 2/(4n) d / 2 T(d/2), z(t) = p(t)/k and U%($) = U%($) /k 2 . Thus, with a given 
Pk,e(p) = t, we first invert the expression to obtain z(t) and then make a substitution into 
Eq. (2.5) to deduce the corresponding smooth RG equation for £/*;($). 

As stated in the Introduction, when a truncated polynomial expansion is used, physical 
quantities such as the critical exponents may depend on M as well as on a. For a sharp 
cutoff, the critical exponent v has been shown to exhibit oscillatory behavior about its 
expected value when £/&($) is expanded about the origin or the /c-dependent minimum 
[7] [12], although the result improves significantly in the latter [9]. In order to determine 
whether these oscillations can be removed or reduced by using a smooth smearing function 
we consider below three parameterizations of Pk,a{p), all of which approach Q(k— p), when 
the appropriate limit of the smoothness parameter a is taken. 
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A. Hyperbolic Tangent 



We first examine 



Pk,e(p) = g 



1 + tanh 



k 2 — p 1 
pke 



The above smearing function satisfies: 



(i) limpfc, e (p) = 0(fc-p), 

e^O 

(ii) lim pk,e(p) = 7i[ 1 + tanh(-oo)] = 0, 

k— >0 Z 

(iii) lim p k ,e(p) = 1 - 

k— >oo 



In this case the propagator is modified as: 



A fc , e (p) 



1 1 

"2 ' 2 



1 + tanh 



Substituting 



dp 



^ 2 + P 2 , 

sech 



2kpe 



kpe 



into Eq. (1.6) then yields: 
0*7 fc ($) 



k- 



J \p k) \ kpe Jp 2 +{1- Pk , 



dk Ae 

With the help of Eq. (2.4) we can rewrite the above expression as 

At) um 



dk 



where 



z(t) = 



z 2 {t) + {l-t)U'*{<S>) 
ea(t) 



~ 2 a 2 (t) 



+ 1 



4 2 
a(t) = tanh _1 (2t - 1). 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



One of our main goals here is to examine how the flow of the theory is influenced by tuning 
e. As e — > 0, the expected sharp-cutoff Wegner-Houghton equation [6] 



k 



dU k {<S>) S d k d 



dk 



dt 



S d k d 



2 j l+ii-tmw 



In 



fc 2 + um 
k 2 



(2.13) 



is recovered. As e is increased the cutoff becomes smoother, i.e., the peak in dpk, e (p)/dp at 
p = k spreads. At about e xs 2 this tendency is reversed and the peak becomes taller again. 
The location of the peak stays roughly the same until e ~ 1 and is shifted towards smaller 
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values as e is further increased. It finally moves to p = k/e for large e. The decrease of 
the cutoff when e — > oo eliminates the interactions and the theory flows into the trivial 
Gaussian fixed point. We shall comment more on this point later. 
Notice that the mean approach described in Eq. (2.1) gives 

k dk 4eJ dpP {p + k) SeCh \ hps )p*+[l-PK,e(p)]um 

Sd r , d-i M(h , g£(g) S d k d ( V»i$) \ 

(2-14) 

which differs from the Wegner-Houghton equation in its absence of the characteristic log- 
arithmic functonal structure. In fact, the RG equation is highly sensitive to the k depen- 
dence in Pk,a(p)- Nevertheless, we shall see that the difference between the two equations 
does not affect the critical properties significantly in the next section. 

Next, we turn the nonlinear flow of £/&(<&) into a set of coupled ordinary differential 
equations by making an expansion in power series of <f>: 



- (2Q! 



followed by a truncation at some order £ = M. The resulting running equations for the 
first two terms read: 



du 2 

- -2 

~dk = h ' k 



(2.16) 



where e = 4 — d, 



J n <jil) = -^L f dt ^ (t)d+2(n+1) 2+ , (2.17) 

and z(t) is given by Eq. (2.12). In the sharp cutoff limit where z(t) — > 1 and 

Jn(fil) = TTTT^-T \~ 2 ^ + 3ll\ (2-18) 



2(1 + n) nKl + fil) 1 * ' 



we recover 



u d »l _ or, 2 Sd Xk 



2 ' 



dk ™ 2 1 + p, 2 

k^ = -e\ k + ^f 3X1 ^ 



dk " 2 l(l + ^) 2 1 + 



(2.19) 
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This can be compared with the result obtained using the mean approach: 



/-. fyfc _ 2 „2 Sd 



2 ' 



-J6) . ( 2 - 2 °) 



(l + «) 3 (l + «) 2 



A comparison of the two flow equations shows that in order to recover the standard per- 
turbative RG coefficient functions in the UV limit ^<Cl, the usual convention 6q = 1/2 
must be imposed. However, this choice also implies a deviation of scaling for a massive 
theory in the IR regime where ji 2 3> 1. 

It is easy to understand that the two flows agree in the UV scaling regime only. In 
fact, the non-uniform convergence in Eq. (2.2) leads to an "incorrect" use of the integrand 
of the evolution equation Eq. (1.6) at p ~ k when the limit e — > is taken before the 
integration. Since the integrand is independent of the field $ for k 2 3> U%($>) the mistake 
is unimportant. But the naive evolution equation displays wrong U{!(<&) dependence when 
k 2 w U%($). 



B. Exponential Function 

An alternative class of smearing functions is the exponential function 

PkM = e-°(*/*>\ (2.21) 
where a and b are constants. This smearing function satisfies the following conditions: 



(2.22) 



In order to facilitate comparison with the other smearing functions used here, we also 
require Pk,b(p = k) = 1/2, which in turn implies a = In 2, or pk,b(p) = 2~( p / fc ) . 
The corresponding RG equation reads 

= at r dy 2 r^-,, - v (2.23) 

Alternatively, one can also cast Eq. (2.23) into the same form as Eq. (2.11), but with z(t) 
given by 

^) = (-^r. (2.24) 



(i) 


lim p k ,b(p) 

b — »oo 


= e(k- P ), 


(ii) 


limp fej6 (p) = 


= 0, 


(iii) 


lim pkM 

k — >CXD 


= i. 
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In the asymptotic limit where U'^{<$>) <C 1, the RG flow can be approximated as 



dk 

where 



k dU k S d fr 



J J" TT" 2 77" 3 ' 



(2.25) 



Cl = a-( d - 2 )/ b r( 6 + ^~ 2 ), 

, , , (2.26) 

C2 = _[i_2-w-4)/t] -(d-4)/b r( 6 + *~ 4 ). 

The sharp limit, on the other hand, gives ci = 1 and C2 = —1/2 independent of dimen- 
sionality d. Thus, c\ and C2 provide a measure of the deviation from the sharp-cutoff UV 
scaling. We also note that in order to avoid UV singularities in the integrations, b must 
be chosen as to avoid singularities in the Gamma functions above. As b is decreased the 
cutoff becomes smoother. The location of the peak in dpk,b(p)/dp is stable for large b and 
starts to shift at b ~ 3. For b < 3, the shape of Pk,b(p) spreads in such a manner that one 
can no longer associate a well-defined cutoff value. 

C. Power-law Function 

Lastly, we consider a power-law smearing function: 

P*.m(p) = 1+( £ /Jfe)m > ( 2 - 27 ) 

which satisfies 

(i) lim p k , m (p) = Q(k -p), [6 = 1/2 when p = k], 

m — >oo 

(ii) lim Pk,m(p) = 0, (2.28) 

(iii) lim pk,m{p) = 1- 

fc^OO 

The RG equation in this case reads 



dk 2 Jo {l + y m )[y 2 {l + y m ) + y m 

S d k d f 1 dt z d (t) UZ(Q) 



(2.29) 



2 J z*{t) + {l-t)UZW 
where y = p/k and 

f \ n 1/m 

z(t) = (- " l) • (2-30) 
The sharp limit, i.e., z — > 1, may also be obtained by the substitution: 

>2 G(—^,p)^6(l-y) I dtG(t,p), m^oo, (2.31) 
(l + y m ) L ~ry Jo 
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which results from 



lim 

m^oo 



my 



TO — 1 



(1 + 2/ 



m\2 



= S(l-y). 



(2.32) 



The corresponding flow equations for n\ and Xk have the same functional forms as those 
in Eq. (2.20) but with 



S f°° 

Jn(fik) = f m J o dz 



^ _|_ ^m^n gd+2n+m+l 

[z 2 (l + z m )+z m fll 



_ 91 2+n " 



(2.33) 



On the other hand, if one adopts the mean approach and takes the limit m — > oo, the RG 
equation then becomes 



k 



dU k {$) S d k< 



dk 



dy 6(1 -y) 



y d u' k 'm 



y2 + e(y-l)C^(*) 



S d k d ZTgg 
2 k 2 + 6 U"(Q) 



, (2.34) 



which again coincides with Eq. (2.14). The peak in dpk, m (p)/dp is stable at p = k for 
m > 5 but shifts towards lower momenta for m < 5, and the shape of the smearing function 
makes it difficult to identify a clear cutoff value. 

Unlike the previous two parameterizations where the integro-differential flow equations 
must be solved numerically, for certain values of d and m analytic evaluation of the integral 
on the right-hand-side of Eq. (2.29) is possible. For example, taking d = m = 4 we have 



k 



dUk 
dk 



k 4 



16tt 2 



a l _ 2 tan -i r k \ \ 

\M 4 - uf I V 4fc4 -^ 



(2.35) 



And for d = 3, we have 



ok 



47T 

k 3 



4tt y/2 + UZ($)/k 2 



-2+J4+2U' t !(<S>)/k< 



(m = 2), 



(m = 4). 



(2.36) 

In this differential form it is possible to apply the standard techniques of differential equa- 
tions to solve for the RG flow of the theory. 

In the asymptotic limit where <C 1, Eq. (2.29) becomes 



k 



dU k 
dk 



S,k d 



U 



//3 



k G 



(2.37) 



where 



m 



ci = 
c 2 = 



d + m-2 

m 



r( 



m 



in 



1, 



2(d + 2m-4) 



r ^m-d + 4^3m-d + 4^ 



m 



1 

2" 



(2.38) 
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Before closing this section, we remark that the RG techniques discussed so far can 
be readily extended to finite-temperature systems. For the scalar A$ 4 theory defined on 
S 1 x R d in the imaginary-time formalism, with the radius of 5" 1 being given by /?, the 
inverse temperature, one obtains 



k 



dk 



S d k d 

w 

S d k d 

P 



P^k 2 + U'^ k ($) + 2 In [l - e -N k ' 2 + u l^ 



In sinh 



(2.39) 



using a sharp cutoff [16]. A smooth smearing function, on the other hand, yields an 
integro-differential equation: 



II" 

U /3,k 



9k 4J P { ^ 1 ^ 2 +{^-PKe)Ul k 



coth(^ ) 



S d k d 



t(z=oo) 



Jt(z= 



dt 



*<*=o) Jz*(t) + (l-t)Ug tk 



(2.40) 

where f3 = f3k. The two equations coincide in the sharp cutoff limit. In the high- 
temperature regime where ft — > 0, we have 



, dU 0k (&) S d k d r^=™) 
k l ^- L « -Ar- I dt z d (t) 



dk 



t(z=0) 



(2.41) 

Apart from the factor of (3 , the first term in the bracket corresponds to the zero- 
temperature (i-dimensional case, in accord with the expectation of dimensional reduc- 
tion. In this limit, j3 can be scaled away with a redefinition of Up^i®) — > /3~ 1 Uk(&) and 
$ [16] [19]. 



III. FIXED-POINT SOLUTIONS AND NUMERICAL RESULTS 

As seen in the last section, when the smearing function is smooth the flow of the 
theory is generally characterized by an integro-differential equation, as opposed to the 
differential equation obtained with the sharp cutoff or special cases of a smooth cutoff. In 
dimensionless form, the flow equation reads: 



-f 

JO 



dt 



(3.1) 
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where U k ($) = ( 2 k- d U k ($), $ = C*:~ (d ~ 2)/2 $, C = VW^ = ^/(4n) d / 2 T(d/2), and the 
momentum scale z(t) is given by Eqs. (2.12), (2.24) or (2.30) for the three cases considered. 
That Sd can be absorbed by a redefinition of $ is an indication of universality, i.e., critical 
exponents are independent of Sd, though the matrix elements in the linearized RG matrix 
are. In the above, the anomalous dimension r\ has been set to zero since we are only 
considering the flow of the potential. 

To characterize the critical behavior of the theory, our first task is to identify all the 
fixed points and then linearize RG about a particular fixed point. It is well known that 
for the one-component scalar theory in d = 3, there are only two fixed points: one trivial 
Gaussian and one Wilson-Fisher, and the critical behavior is dominated by the latter. No 
other continuum limit is known to exist. This applies to both the sharp and the smooth 
cutoffs. However, for a general smooth function Pk,a(p), the location of the Wilson-Fisher 
fixed point which will now depend on a; at a = 0, however, it coincides with the Gaussian 
one. 

At the fixed point(s), the theory exhibits scale invariance, i.e. dkU£ = 0, and the RG 
equation for the fixed point potential becomes (dropping the subscript k): 



1 f 1 

■-(d-2)W*'($) + dU*{$) = - \ 

2 Jo 



dt ^ U * ^ 



" -In 



z*(t) + (l-t)U*"(*) ( 3. 2) 
1 + &*"($) 



where the notation s. stands for the sharp cutoff limit. Analytically, the sharp limit has 
the following approximate non-trivial solution in the large $ limit for d = 3: 

U*($) = - - In $ - - - - ln(3(L4) + 0($" 6 ), (3.3) 

v 1 3 9 3 v 1 15(L4$ 4 V ; V 1 

where A is an arbitrary positive constant [12]. The RG flow about Eq. (3.3) can be 
linearized by writing U k (&) = U*($>) + v(§)e~^ lnk , where v($>) obeys 

--{d - 2)$v'{$) - (A - <£)v($) = - [ dt z d (tW'($) 

2 ^ Jo [z*{t) + (l-t)U*"{*)] 2 (34) 

s. v"(Q) 

— ► ~ i + u*"(®y 

We comment that the location of the fixed point as well as U *($>), being non-universal, 
naturally all depend on the choice of the smearing function Pk(p)- In the one extreme where 
Pk,a(p) = ®(k —p), linearlizing RG and retaining only the relevant operators yields 

(^,A«) = (0,0), (--i-,-^), (3.5) 

where B = 1/16/t 2 . However, in the other extreme where the smearing function is taken to 
be a /c-independent constant, Pk,a(p) = c, the right-hand-side of Eq. (3.1) then vanishes, 
and the RG flow simplifies to 

dUk^) = Q ^ -l(d-2)®U*'($) + dU*($) = 0. (3.6) 
dk 2 
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The equation can be exactly solved to give U*($) = A^ 2d ^ d ~ 2 \ For e = 1 or d = 3, 
U* ~ $ 6 . Thus, with vanishing contributions from quantum fluctuations, the only fixed 
point is Gaussian (/i 2 *, A*) = (0,0). However, this does not mean that the theory is free, 
because we still have U*(&) ^ 0. 

For the general case of a smooth cutoff function, however, it is rather difficult to obtain 
a non-truncated solution for Eq. (3.2), and we circumvent the problem with polynomial 
truncation at order (j) 2M . The remaining task is to solve a system of M integro-differential 
equations for the coupling constants g^ 2 ^*, £ = 1,---,M. We first determine the fixed- 
point solution analytically in the sharp cutoff limit, so that we can then use this as an 
initial guess in our numerical root-finding subroutine for the smooth case. The system of 
integro-differential equations is solved by first choosing a large value for the smoothness 
parameter a 6, or m), thereby making the smearing function rather sharp. The 

smoothness of the cutoff function is increased in small steps by decreasing the parameters, 
and the solution of the previous step is used as the guess for finding the solution at the 
current step. In this manner we are able to track the solution associated with a physical 
fixed point from the sharp cutoff limit to an arbitrary smoothness. 

With the strategies outlined above, we now present the results for all three smearing 
functions considered in Sec. II. 
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Figure 1. Critical exponent v as a function of 1/e using the hyperbolic tangent smearing function. 
Results for four different levels of polynomial truncation are shown. 
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In Fig. 1 we plot the dependence of v in d = 3 as a function of the (inverse) smoothness 
parameter e -1 at four different levels of truncation of [/*(<&). From the figure, we see that 
v varies by 2-5% over the range shown. The variation would become even larger had we 
taken e~ x — > 0, where the smearing function becomes very flat and there is practically no 
blocking. The exponents in this limit are those obtained in the mean-field approximation, 
e.g., v = 0.5. It is also interesting to note that at each order in M there is a dip, or 
an extremum whose position changes only slightly with M. The non-monotonic behavior 
shows that when going to the sharp limit by gradually increasing e -1 , v first deviates from 
the sharp result and then begins to converge when pk,e(p) becomes sufficiently sharp. 

When an exponential smearing function is employed instead, we observe a similar 
behavior for u, as depicted in Figure 2. Once more, the value varies by 2-5% and an 
extremum is found for each M. The qualitative feature remains the same for a power-law 
smearing function, as demonstrated in Figure 3. However, in this case the variation of v 
is only about 2-3%. 
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Figure 2. Critical exponent v as a function of b for the exponential smearing function. Results for 
four different levels of polynomial truncation are shown. 
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Figure 3. Critical exponent v as a function of m for the power law smearing function. Results for 
four different levels of polynomial truncation are shown. 

IV. SCHEME DEPENDENCE 



Even though non-truncated solutions for the fixed-point potential £/"*(#) can be ob- 
tained in special cases [12], this is not possible in general. Any approximation inevitably 
brings in scheme-dependent effects. Our results are seen to depend both on the order of 
truncation M as well as the smoothness parameter a. We discuss below the implications 
of such dependence. 

A. Truncation Dependence 

At order M we obtain M fixed-point solutions parameterized by the critical coupling 
constants (g^* ■■■ g^ 2M '*), each having its own eigenvectors and eigenvalues which in 
turn can be used for calculating the critical exponents. However, only one of the solu- 
tions corresponds to the Wilson-Fisher fixed point which is characterized by one relevant 
operator, with the rest being irrelevant. The trivial Gaussian solution with all g( 2i ^* = 
is also obtained but is of no interest to us. Computational artifacts are induced with a 
polynomial truncation of the blocked potential Uk(&) [9] [12] [16]. For example, as M is 
increased, numerous unphysical fixed points are also generated. 

Thus, to isolate the physically meaningful solution from the unphysical ones in the 
smooth parameterization, what we have done numerically was to first analyze the sharp 
cutoff limit of Eq. (3.2) and identify the solution associated with the Wilson-Fisher fixed 
point. The salient feature of this solution is that it is stable against an increase in M. In 
other words, the value of g^* is not sensitive to whether higher order coupling constants, 
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e.g., g( 2 ^+ 2 )*, are included or not. On the other hand, unphysical solutions are sensitive to 
the higher order equations and their eigenvalues will fluctuate in an unpredictable manner 
[8]- 

Another prominent artifact of the polynomial truncation is the oscillation of the 
critical exponents with M when a sharp cutoff is used. The values of the critical cou- 
pling constants also alternate in sign. The oscillations are due to the nonanalyticity 
of the critical blocked potential [/*($) and can be understood from the fact that an 
expansion of the right-hand-side of the sharp cutoff flow equation ln(l + U* ($)) = 
J2iLi(— 1)^ +I )(f7* (&))*/£ also yields a power series of $ with alternating sign. Even 
though the convergence has been shown to improve significantly by expanding 
around its moving minimum, there still remains residual scheme dependence [9]. In addi- 
tion, an expansion about the minimum fails in the broken phase of O(N) as stated in the 
Introduction. Therefore, it is desirable to have a numerical technique which allows for an 
efficient determination of the critical properties, and remove or reduce the aforementioned 
spurious effects. 

Another subtle issue remains in this approximation is: why can we not reproduce the 
exact critical exponent(s) by keeping only the renormalizable coupling constants in our 
computation and in the identification of the fixed point? The main reason lies in that fact 
the fixed point of the truncated solution is applicable only up to the operators which have 
been neglected. At the point which we call a "fixed point," the would-be ignored set of 
irrelevant operators does not vanish but continues to evolve. On the other hand, in the 
exact RG approach, the fixed-point condition implies a complicated cancellation between 
different operators. When all the operators are present, the cancellation is complete and 
we should be able to eliminate the pfc )(T (p)-dependence in the scaling laws around the fixed 
point as well as in the expression for the critical exponents. 

B. Smoothness Dependence 

In the alternative Wilson or Polchinski RG approach, one can show that when the 
effective action is expanded in terms of derivatives to order n, only 2n parameters are 
required to absorb the scheme dependence [15]. Therefore, at the leading order of Uk(&) 
with n = 0, the critical exponents must be scheme-independent, i.e., the same results are 
obtained regardless of whether the smearing function is sharp, exponential or power-like. 
In fact, Pk,a(p) can be conveniently absorbed by a suitable redefinition of the field variable 
[20]. Scheme depedence appears, however, when the wavefunction renormalization Zfc(<&) 
is taken into consideration at the next order. 

On the other hand, Eq. (3.1) derived with LPA does not possess this scaling property, 
and the functional form of the flow equation is explicitly dependent of Pk,a(p)- The depen- 
dence, nevertheless, is expected to be negligible in the calculation of physical quantities 
such as the critical exponents. This universality hypothesis, however, must be substan- 
tiated by solving the integro-differential equation Eq. (3.1) exactly without polynomial 
expansion. 

From Figures 1, 2 and 3, the critical exponent v clearly depends on the smoothness 
parameter a (e -1 , b and m). This comes as a consequence of polynomial truncation. That 
is, because of the polynomial truncation scheme employed, physical results now depend 
both on the level of truncation M as well as on a. 
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How can such dependence be reconciled with universality? The key to the problem 
here again lies in the role played by the irrelevant operators in the approximate solution. 
If we solve the RG equation for the full theory and come down to the IR regime, the 
dependence on the initial condition for the irrelevant operators must die out according 
to the universality hypothesis. The explicit verification of this scenario by reaching the 
IR fixed point in principle requires an infinite number of iterations. This is not what 
was followed in this work. Instead we move close to the UV fixed point and a linearized 
RG prescription is utilized to deduce the critical exponents and the corresponding scaling 
laws. Using the derivative expansion and the polynomial truncation schemes gives the 
advantage of a readily accessible (approximate) fixed-point solution; however, we are no 
longer certain what the omitted irrelevant operators do. If they are not scale invariant at 
the fixed point - the likely result of the truncation, then their variation amounts to the 
change of the overall scale factor of the exact solution of the theory. This is because the 
irrelevant coupling constant set of the theory influences an overall scale factor. In order to 
minimize this scale dependence the truncation scheme must be "improved" . 

The comparison with the exact RG offers a direction of improvement in our approx- 
imation: to minimize the number and the strength of the neglected irrelevant operators. 
As far as the pA; )(T (p)-dependence is concerned, the momentum dependent, i.e. higher order 
derivative terms come into consideration which are quadratic in the field Therefore, we 
aim to minimize the number of such terms generated in the RG flow. The cutoff-function 
Pk,e{p) is strongly momentum dependent in an interval Ap = ek. According to the un- 
certanity principle this yields non-local interactions on the length scale Ax = 1/ek. In 
order to eliminate such a higher order derivative effect from the theory we must look for a 
cutoff which is as smooth as possible in the momentum space. This would yield a better 
convergence when the irrelevant derivative operators are truncated, in accord with the 
spirit of LPA. 

However, there are also higer-order field operators and we should keep M as low as 
possible. We emphasize that the purpose of employing a truncated polynomial expansion 
of is to devise a simple algorithm which allows us to reproduce the nonperturbative 

features, particularly the fixed-point structure embedded in LPA, or even the blocked 
action if?fc[<&], without having to solve the full flow equation for [/*.($). Thus, to go to a 
large M not only is computationally laborious but also in conflict with the original goal. 

In Fig. 4 the truncation dependence of v obtained using the exponential smearing 
function given in Eq. (2.21) is depicted. The sharp cutoff results are also included for 
comparative purpose. We notice the general trend of a more pronounced oscillatory be- 
havior as b is increased toward the sharp limit. This again show that the sharp cutoff is 
not the best candidate if we look for the smallest set of operators which yields the most 
rapid convergence for the physical observables. Below we demonstrate how this can be 
achieved with an optimization scheme. 
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Figure 4. Critical exponent v as a function of M using the exponential smearing function. Results 
for several values of b are shown for comparison. 



C. Optimization 



For a sharp cutoff, with a clear boundary between the high and the low modes, the 
IR cutoff is determined unambiguously as k. However, we wish to choose the cutoff to be 
as smooth as possible in order to minimize the generation of non-local, higher-derivative 
irrelevant interactions, as well as the order of truncation M. A general smooth smearing 
function pk,a(p) results in a shift in the actual value of the IR cutoff, or the location of 
the peak of —dpk, a (p)/dp. When Pk,a(p) is too smooth, we no longer have a well-defined 
cutoff, i.e., the relation between the parameter k and the momenta of the modes appearing 
in the evolution equation. 

As we have shown before, with the truncation scheme employed, the fixed-point solu- 
tions will vary with M and the would-be irrelevant operators will continue to contribute 
and evolve around the these approximate solutions. Thus, we would not expect a priori 
an accurate result for the critical exponents in this case. Nevertheless, we notice that for 
each of the three smearing functions discussed in Sec. II, there is an "optimal" value of 
a for which v converges most rapidly. For example, the optimal value for the exponential 
function is found to be b = 3, which incidentally coincides with the value of the extremum 
shown in Figure 2. When all three smearing functions are superimposed, it can be seen 
that they have approximately the same shape or smoothness, as illustrated in Fig. 5. That 
is, the most optimal smoothness does not depend on the detailed form of the smearing, 
i.e., hyperbolic tangent, exponential or power-law. 
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Figure 5. Comparison of the three optimized smearing functions. 

The truncation dependence of v for the three optimized smearing functions is ilus- 
trated in Fig. 6. Although there remains small oscillations, we see a dramatic improvement 
in the convergence. In fact, the variation beyond M = 7 is only 1-2%. Our optimized value 
gives v = 0.65(5) which is closer to the world's best value 0.631(2) [21] than 0.68(9) ob- 
tained by Aoki et. al. using an expansion around the moving minimum. We also remark 
that M, the number of operators involved in our calculation, is considerably smaller for a 
reasonably accurate estimate of v [9]. However, by truncating the potential at M = 12, 
we have not been able to observe the oscillation of v with M at a four-fold periodicity; an 
unambiguous observation of such behavior generally would require M > 20 [9] [12]. 
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Figure 6. Critical exponent v as a function of the level of polynomial truncation for all three opti- 
mized smearing functions and the sharp cutoff. 

Why is a related to the convergence of the series expansion? Physically, the param- 
eter dictates the manner in which the quantum fluctuations modes are integrated over. 
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When a is very small and the corresponding smearing function is too smooth, certain fast- 
fluctuating modes are only damped but not integrated over due to the lack of a well-defined 
boundary. In addition, some slowly-varying degrees of freedom which should be kept are 
integrated. As a result of the "mistreatment" of both the fast and the slow modes, the RG 
trajectory is greatly distorted, and the would-be irrelevant higher-order operators continue 
to evolve near the critical point obtained at order M. 

The complicated interplay between the irrelevant operators can also be seen from the 
non-mo notonic dependence of v on a observed in Figures 1 to 3. An extremum is found 
at cr c for all orders of M. While for a c is a local maximum for M = 2, it becomes a local 
minimum for M > 3. The flip is clearly due to the inclusion of the irrelevant (marginal by 
power counting) $ 6 operator. 

In general, as a increases and more fast modes are being included in the loop inte- 
grations, more cancellations between the effects of the irrelevant operators take place and 
the theory moves closer to the true Wilson-Fisher fixed point found in the exact approach. 
Nevertheless, when a becomes too large and the resulting smearing function is too sharp, 
non-local effects begins to set in and the theory drifts away from the true RG trajectory. 
At cr = oo where the smearing function becomes a step function Q(k —p), non-local effects 
become maximal and LPA is no longer adequate. One must then include derivative oper- 
ators to all orders in the blocked action Sk[&] in order to arrive at a scheme-independent 
result. 

Thus, we see that a monitors the manner in which the irrelevant operators contribute 
to the flow of the theory. When the most optimal smoothness is reached, the maximum 
cancellation between these operators takes place, leading to the fastest convergence in the 
polynomial truncation. 

In Sec. II we have also discussed an alternative mean approach, which is based on 
taking the sharp cutoff limit followed by a substitution of 6$ = 1/2 to avoid ambiguity, 
and yields Eq. (2.14) as the RG flow. We solve the equation and depict the result for v in 
Figure 7. 

Contrary to the sharp approach which results in an oscillatory behavior, v is seen to 
converge rapidly. Such a change in the convergence is intimately related to the structure 
of the untruncated fixed-point potential £/*(<&) [5] which now obeys 
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Figure 7. Critical exponent v as a function of the level of polynomial truncation obtained with the 
mean approach and the sharp cutoff. 



Even though Eq. (4.1) still cannot be solved analytically, its improvement in con- 
vergence can be shown by examining the right-hand-side of the equation representing 
the quantum corrections. By making a Taylor expansion of the expression with respect to 
we see that it has a radius of convergence 2/Uj{ which is twice that of ln[l + £/£'($)] 
in Eq. (3.2). Since a general smooth cutoff always has a greater convergence radius, the 
resulting critical exponents will converge more rapidly compared with the sharp in the 
polynomial truncation scheme. To describe the amplitude and the periodicity of the oscil- 
lation more quantitatively, however, would require a detailed knowledge of the singularities 
of the untruncated potential £/*($). 

V. SUMMARY AND DISCUSSIONS 

In the present work we have demonstrated how one can improve the convergence of 
critical exponents calculated using polynomial truncations of the RG flow equation for the 
blocked potential Uk(&)- Since in the sharp limit, v oscillates with M, the order of trunca- 
tion, we examined three parameterizations of a smooth cutoff function, hyperbolic tangent, 
exponential, and power-law, and made an attempt to eliminate as much as possible such 
unphysical artifacts. We find that there exists an optimal smoothness value which gives 
the most rapid convergence as a function of M. This is due to the maximal cancellation 
of the effects generated by the irrelevant operators in the blocked action. Our optimal 
smearing functions yield v = 0.65(5) for M > 7 with a variation of 1-2% between the three 
cases. Had a been too large, we would have to incorporate the higher-order derivative 
operators to account for the non-local effects generated in the course of RG evolution. 

We have also learned from comparing the results for the sharp limit and the mean 
approach (depicted in Figure 7) that the convergence behavior of the polynomial truncation 
scheme is intimately related to the non-truncated solution of (<!>). The mean approach 
results in a larger radius of convergence, and hence a more rapid convergence. It remains 
to see if the same "trick" can be utilized for more complicated systems as well. 

In light of the success of our optimized RG prescription, we can readily extend the 
formalism to address other issues such as the O(N) models, the spinodal instability [22], 
gauge theories or chiral symmetry breaking. It would also be interesting to compare our 
usual momentum-cutoff approch with an alternative internal-space RG, which is a func- 
tional generalization of the Callan-Symanzyk equation, based on an infinitesimal variation 
of an internal-space parameter such as the mass scale [23]. 
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